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1. Introduction 

The equations of the relativistic magnetohydrodynamics (RMHD) of a magnetoactive 
plasma in a gravitational field were formulated in [lJJ using the equality requirements 
for the dynamic velocities of the plasma and the electromagnetic fielq§|. These 
equations were obtained on the basis of the Einstein and Maxwell equations. A 
remarkable class of exact solutions of these RMHD equations was also found. It 
explains the motion of a magnetoactive locally isotropic plasma in the field of a 
plane gravitational wave (PGW). This class was called gravimagnetic shock waves 
(GMSW). It describes essentially nonlinear processes which do not exist in the linear 
approximation of magnetohydrodynamics and essentially relativistic processes in terms 
of predominance of the massless electromagnetic component in the magnetoactive 
plasma. It was shown in [2] that the GMSW in pulsar magnetospheres may be the 
highly effective detectors of gravitational waves from neutron stars. Particularly, giant 
pulses which sporadically appear in the radiation of some pulsars may be observational 
results of energy transfer from a gravitational wave to GMSW. Estimations made 
in [2]- [4] make it possible to connect giant pulses in radiation of the pulsar 
B0531+21 with gravitational radiation in the basic mode of oscillations from this 
pulsar. In fact, at present it is rather difficult to speak of identification of giant 
pulses as an electromagnetic display of gravimagnetic shock wave evolution in the 
pulsar magnetosphere and to unambiguously connect these pulses with the pulsar's 
gravitational radiation. Nevertheless, the idea of analyzing the effect of gravitational 
waves from a compact astrophysical object on its own electromagnetic radiation is 
highly productive for solving the problems of gravitational waves detection. In fact, 
the main difficulties of gravitational waves detection in the Earth environment are: 

(i) An extremely small amplitude of gravitational waves on Earth (h ~ 10" 19 ) due 
to significant distances from relativistic astrophysical objects. 

(ii) A sporadic nature of events leading to radiation of gravitational waves inside 
relativistic astrophysical objects with sufficient power. This does not allow one 
to unambiguously connect a received signal with a fact of gravitation radiation 
detection. 

(iii) Impossibility to construct relativistic detectors with anomalous, highly effective 
parameters for gravitational wave detection in conditions terrestrial laboratory 
(superstrong magnetic fields, a highly anisotropic working body of the detector, 
low level of background noise etc.). 

It is possible to avoid these problems if one could transfer a detector directly to a 
close neighborhood of a relativistic astrophysical object. In this case, one always 
has a ready electromagnetic signal and there is no need to convert it to other forms, 
which allows for conducting correlation analysis. If the working body of the detector 
is the magnetosphere of a relativistic astrophysical object, the optimal parameters 
for gravitational wave detection are achieved automatically: super-strong magnetic 
fields, an ultrarelativistic equation of state, highly anisotropy etc. The fundamental 
importance of the GMSW for gravitational theory, as a direct conversion effect of 
gravitational wave energy into electromagnetic energy, leads to a necessity of a more 
detailed and comprehensive researches. In |5,, a strict foundation of the GMSW 

\ Before 2000 Yu.G. Ignatyev wrote his name as Yu.G. Ignat'ev. 

§ This requirement is completely equivalent to the condition of plasma infinite conductivity, see 
Ref.[l]. 
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hydrodynamic theory was formulated on the basis of the relativistic kinetic theory. 
As has been shown in PQ-[I], a GMSW is realized in an essentially collisionless 
nonequilibrium plasma in anomalously strong magnetic fields. The isotropy of a local 
plasma electron distribution essentially is violated under such conditions due to strong 
bremsstrahlung. Therefore, the anisotropy factor of a magnetoactive plasma is highly 
essential for the efficiency of the GMSW formation mechanism. A hydrodynamic 
model of GMSW in an anisotropic plasma was constructed in [6] for a specified relation 
between parallel and perpendicular components of the plasma pressure. The analysis 
in [B] was based on the general RMHD equations. Particularly, an elementary linear 
relation was considered in j6]. This study has revealed a strong dependence of the 
GMSW effect on the plasma degree of anisotropy. That fact has led to a necessity of 
constructing a dynamic model for the motion of an anisotropic magnetoactive plasma 
in the field of gravitational radiation. 

Further, in [7] a detailed numerical model of GMSW has been carried out at 
various parameters of the anisotropic magnetoactive plasma in the computer algebra 
system Mathematica. This study was based on the numerical solution of the nonlinear 
energybalance equation by means of special methods of numerical integration. The 
results received in [7] have confirmed earlier made analytical estimations of the 
magnetoactive plasmas behavior in a field of strong gravitational wave and have 
defined more exactly some characteristics of GMSW. 

However, at all variety of models of plasma only the case of monopolarized 
gravitational wave with polarization e + investigated in all quoted papers. It has been 
thus shown that the case of monopolarized gravitational wave with polarization e x is 
reduced to the case of e + polarization state at the coordinate and physical quantities 
transformations. However, the case of simultaneously existence of two polarizations 
of a gravitational wave was not investigated. In this paper we consider such a case. 
Thus it was possible to reduce some additional conditions, which in [1] have defined 
the structure of potential electromagnetic field in a magnetized plasma, and thereby 
give a more general meaning to results obtained in preceding papers. In this paper 
adopted a system of units where (c = G = %= 1). 



2. Self-consistent RMHD equations in a gravitational field 

2.1. Frozen-in condition of magnetic field in plasma 

In pQ under the assumption of equality of dynamic timelike velocity v l of a plasma 
and electromagnetic ficlcUj]: 

v f 

TijV J =e p Vi] Ti 3 v 3 = EfVi] (v,v) = l (I) 

on the basis of conservation of the total energy-momentum tensor of a plasma and 
electromagnetic field, 

v ■ f ■■ 

T*J = (3) 

full self-consistent system of relativistic magnetohydrodynamic equations for 
magnetized plasma in arbitrary gravitational field has been obtained. It describes the 

| The index "p" refers to the plasma, the index "/" to the field, a comma denote covariant derivatives. 
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motion of a relativistic plasma and an electromagnetic field in the given gravitational 
field. 

In particular, it was shown that at positivity of the first invariant of an 
electromagnetic field: 

Iiwi = FijF^ = 2H 2 > (4) 

and equality to zero of the second invariant: 

Inv 2 =Fij F ij = (5) 

necessary and sufficient condition for the solvability of equations ([T]) is the frozen- 
in condition of magnetic field in plasma, i.e., equality to zero of the accompanying 
intensity of the electric field Ei : 

Ei = F jiV j = 0. (6) 

In formulas (j4j-([6]) and further: Fij - antisymmetric Maxwell tensor, Fij - dual 
Maxwell tensor. 

Fij= \m 3 kiF k \ (7) 

where rjkilm ~ covariantly constant discriminant tensor (Levi-Civita's tensor) [5]. 

At the frozen-in condition ([5]) fulfilment, the condition of dynamic velocities 
equality (JXJ) is always satisfied, regardless of the conditions (@|-([5]). On the basis 
of rigorous kinetic model of plasma one of the authors has shown that the frozen-in 
condition is a consequence of the drift approximation, i.e., the smallness of Larmor 
length for electrons A e = c/w c in comparison with the characteristic inhomogeneity 
scale, r: 

A= «1; lo c = — , (8) 

ruj e m e 

where w c is Larmor frequency for electrons. 

2.2. Self- consistent equations of magnetohydrodynamics 

The complete system of self-consistent RMHD equations for a plasma in a gravitational 
field, obtained in pQ, consists of Maxwell equations of the first group: 

*F l \ k = 0; (9) 
Maxwell equations of the second group: 

F% = -4^J* r (10) 
with spacelike drift current: 

p 

QJpik rpl 

*^dr ~ r? (Jdr,J d r)<0 (H) 

and a conservation law of the total energy-momentum of the system: 

T%+ T% = 0. (12) 

The continuity equation for the drift current must be satisfied in consequence of Eq. 
(HQJ): 

, = 0. (13) 



Exact solutions of relativistic magnetohydrodynamics equations 5 

It should be noted some useful strict consequences of magnetohydrodynamics 
equations: 

Fik Ji = o; (14) 

v l T% = 0; (15) 

//' r:.i, <)• (i6) 

2.3. Maxwell tensor representation by accompanying intensities 

The components of Maxwell tensor is conveniently represented by a pair of spacelike 
vectors of accompanying intensities of electric, Ei ([5]), and magnetic, Hi, fields [5]: 

E i = F ji v j ; Hi=FjiV j , (17) 

so that: 

(E,E) = -E 2 ; (H,H) = -H 2 ; (y,E)=0; (v,H) = 0. (18) 
Then Maxwell tensor and dual to it can be expressed through a pair of spacelike 
vectors of this accompanying intensities [5] : 

F i:i = ViEj - VjEk - i,, j.ir' II 1 : (19) 

Fij= ViHj - Vj H k + r] m v k E l , (20) 



where: 



1 " 1 " - ■ ■ ■ ■ -l2 C 2 



F H Fl3 = o FijF V = E ) " ( ff > H) = H Z -E Z - (21) 



/•'..,/•"' [E,H). (22) 



2 J 2 
1 

2 

The energy-momentum tensor (EMT) of the electromagnetic field 

h = ^ (f\F\ + \5lF lm F lm ^ (23) 

can also be represented by the triplet of vectors v,E,H (see PQ). In the case of 
coincidence plasma's and electromagnetic field's dynamic velocities ([I]) the EMT is 
expressed through a pair of vectors, v,H pQ: 

h = ~ M - 2v l v k )H 2 + 2WH k )] , (24) 

so: 

T=T\ = 0. (25) 
It is easy to verify that the vector v and the spacelike unit vector of the magnetic field 
h - 

h' = ^ (h,h) = -l; (v,h)=0 (26) 

/ ., 

- are actually the eigenvectors of the tensor T ■ 

Tlv k =e H v*; (27) 

T\h k = e H h\ (28) 
where the invariant 

- ff (29) 
is the energy density of the magnetic field. 
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2.4- Energy-momentum tensor of magnetoactive plasma 

The energy-momentum tensor of a relativistic anisotropic magnetoactive plasma in 
gravitational and magnetic fields is (see, for example, [B]): 

T ij = (e+pj>V - P± gV + (p,, - px^W, (30) 

where p_i_,p\\ - plasmas pressure in the directions orthogonal and parallel to the 
magnetic field, respectively. Trace of the energy- momentum tensor (|30D is: 

T=T\ =s- P± -2 P[l > (31) 

and because of the virial theorem (see [10]) it is non- negative: 

P± + 2p|| < e. (32) 

It is easy to verify that the vectors v and h are also eigenvectors of the energy- 
momentum tensor of plasma (see (f2"7f , (|2"51) ) 

T lk v k =ev l - (33) 
T zk h k = - P \ { h\ (34) 

3. Solving RMHD equations in the PGW metric 

3.1. The metric of a plane gravitational wave 

The vacuum PGW metric is (see, for example, |11|): 

ds 2 = 2dudv- L 2 dZ 2 , (35) 

where: 

d£ 2 = cosh27(e 2/3 ((fa; 2 ) 2 -I- e~ 2p {dx 3 ) 2 ) - 2 sinh 2 1 dx 2 dx 3 (36) 

- is a metric of "plane" (x 2 , x 3 ); /3(u), 7(1*) - amplitudes of the polarization e + and e x , 
respectively; u = -^(t~ xl ) is the retarded time, v — -7^ (t + a; 1 ) is the advanced time. 
The amplitudes of PGW are arbitrary functions of the retarded time u, and L(u) is 
a background factor of PGW, which defined by single nontrivial vacuum Einstein's 
equation^: 

L" + L(cosh 2 2 7 /3' 2 + 7 ' 2 ) = 0. (37) 

At inversion of the coordinates in the plane (x 2 ,x 3 ) and transformation of the PGW 
amplitude: 

x 2 =x' 3 ; x 3 =x' 2 ; & = -p\ 7' = 7 (38) 

two-dimensional metric transforms into itself. Under rotations in the plane (x 2 ,x 3 ) 
by the 7r/4 angle: 

x 2 = -±=(x' 2 + x' 3 ); x 3 = -±=(x' 3 + x' 2 ) (39) 

two-dimensional metric is transformed to: 

dS' 2 = (cosh27COsh2/3 + sinh27)(da;' 2 ) 2 + 

(cosh27Cosh2/3-sinh27)(dx' 3 ) 2 + 2cosh27smh2/3efo' 2 dx' 3 . ^° 

% The prime denotes differentiation with respect to the retarded time u. 
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If P — 0, i.e., in the case of PGW with & single polarization g x , we get from (gO) : 

d^ 2 = e 2 ^(dx' 2 ) 2 + e-^(dx' 3 f 

- a PGW metric with a single polarization of e+. 
For a weak gravitational wave: 

|/3(u)|«l; |7(u)|<l; L 2 (u) = 1 + 2 (|/3, 7 |) (41) 

rotation (|39p is equivalent to the transformation of inversion: 

/3'= 7 ; t' = - j 9. (42) 

3.1?. Initial conditions 

Let in the absence of PGW (u < 0): 

/3(w<0)=0; 0'(tt<O)=O; L(u < 0) = 1, (43) 
plasma is homogeneous and at rest: 

v v (u<0) =v u (u<0) = 1/V2; w 2 = i' 3 =0; 

e(it < 0) =e; pj, (u < 0) =J|| ; P± (u < 0) =P_l, (44) 

and homogeneous magnetic field is directed in the (x 1 ^ 2 ) plane: 

o o 
Hi(u < Q) =H cos O ; i? 2 (u < 0) =H sin 17 ; 

H 3 (u <0)-0; Si(u<0)=0, (45) 

where is the angle between the axis Cte 1 (the PGW propagation direction) and the 
magnetic field H. 

As we noted above, the effect of a PGW with polarization e + on a homogeneous 

plasma at the initial conditions (j44"l) -(j45 |) was considered in the quoted papers PQ-[6]. 

Taking into account transformational properties of the metric noted in section 13.11 it 

means that the effect of a monopolarized PGW on a homogeneous plasma has been 

o 

considered earlier when the projection of vector H on the plane of PGWs front is 
parallel to the polarization axis. The case, when this projection coincides with the 
direction x 2 or x 3 , can be reduced to the case of polarization e + or e x using the 
substitution (|38j) . And vice versa: the case of different polarizations e + or e x can be 
reduced to the case with different projections on the direction x 2 or a; 3 under rotation 
in the PGWs front plane by the angle 7r/4 together with rotation of the vector of 
magnetic field intensity. For understanding of a mechanism of strong PGW interaction 
with an anisotropic magnetoactive plasma it is essentially important to consider the 
combined case, when a PGW possesses both polarization states simultaneously, and 
the projection of the vector of magnetic field intensity on the PGWs front plane is 
parallel to the axis of one of them. The initial conditions (|4"4"l) - (|4l)|) correspond to this 
case. 



3.3. Symmetry of the problem 

As is well known, the metric (|35[) permits the group of motions G5, associated with 
three linearly independent in a point Killing vectors: 

r = si ; e = si ; m 

(1) (2) (3) 
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In consequence of the Killing vectors existence in the metric (1351) , all of the geometric 
objects, including the Christoffel symbols, the Riemann tensor, the Ricci tensor 
and, consequently, the energy-momentum tensor of a magnetoactive plasma, are 
automatically conserved at motions along the Killing's directions: 

L 9ij = 0; => L R l3 = 0; L T,j = 0, (47) 
where L is a Lie derivative in the direction of £: 

L Tij = T ijtk £ k + T kj £j + T lk ^ . (48) 

p 

Further we require that tensors of energy-momentum of the plasma Tij and the 
/ 

electromagnetic field Tij inherit the symmetry separately: 

h Tij= 0; (49) 

h Tij= 0; (a = 173). (50) 

Consequences of (f50"|) are: 

L F ti = 0, L Fij= L -ff = 0, L £^ = 0, L ff 4 = 0. (51) 

Consequences of (|4"9")l and (p)Tj) are: 

L £ = 0, L !)' = 0, L pj_ = 0, L pi! = 0. (52) 

Ca «a 



Thus, all observed physical quantities P inherit the symmetry of the metric (|35J): 

LP = 0; (a = T73), (53) 

So, 



i.e., taking into account the explicit form of Killing vectors (|46l) : 

p = p(u); e = s(u); v l = v z (u); (54) 
F ife =F ife (u); Hi = Hi(u); hi=hi{u). (55) 

5.^. Maxwell tensor 

In this section we obtain an expression for the vector potential of the electromagnetic 
field in the metric (|35l) . taking into account the initial conditions (l44l) - (|45l) . This 
method differs from the method used in pQ. It is based only on the first group of 
Maxwell equations and the initial conditions and therefore have greater generality. 
The vector potential conformed with the initial conditions (|4"5|) is: 

A v = A u =^ 2 =0; 

A 3 =H (x 1 sin fl- x 2 cos fl); (u < 0). (56) 
These conditions conform with follow components of Maxwell tensor: 

F 23 (u <Q) = -H cosO; F v3 (u < 0) = —= H smtt; F v2 {u < 0) = 0; 

F u2 (u < 0) = 0; F u3 (u < 0) = --= H sinO; F uv (u < 0) = 0. (57) 

v2 
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As known (see the first group of Maxwell equations ([§]) is equivalent to the 
existence condition of a vector potential. It can be written as: 

1 =d jy /=g F rj = 0. (58) 



Considering (|55|) , we get: 

L 2 F ua = C [a) (= Const); a = {v, 2, 3}, (59) 

setting here and further the following order of the coordinates: 

Coords := [v,u,x 2 ,x 3 ], (60) 

Let us establish a connection between the components of Maxwell tensor with the 
components of tensor dual to it: 

* I * „ I * q 1 

V — —J-^^ 23 ' b — -j-2^v3, t — —-j-2-fv2, 



F v2 = ^F u3 ; p v3 ^^F u2 ; F 23 = ~F m . (61) 

Then the initial conditions (|45"T) give: 

L 2 F w = -F 23 =H cosQ; (62) 

L 2 F 112 = F v3 = 4= H sin fi; (63) 
v2 

L 2 F u3 = -F v2 = 0. (64) 
Thus, the second invariant of the electromagnetic field is equal to: 

Inv 2 = F lk F lk = -^(F v3 F u2 ~ F 23 F UV ), (65) 

So, taking into account (I53")) . the equality to zero of the second invariant of an 

electromagnetic field ([5]) is reduced to the relation: 

L 2 F v3 = F u2 = -V2F UV cot ft. (66) 

As it is known (see for example [TO] ) , the first group of Maxwell equations is equivalent 
to the existence condition of a vector potential Ai : 

F lk = diA k - d k Ai. (67) 

Let us notice that as opposed to Maxwell tensor, the components of the vector 
potential Ai can depend on the variables v,x 2 ,x 3 . We write down expressions (|62|) - 
(|66l) relative to the vector potential Ai using definition of Maxwell tensor (|571) : 

o 



d 3 A 2 


-d 2 A 3 


d v A 3 


-d 3 A v 


d v A 2 


~d 2 A v 


d u A 2 


-d 2 A u 



^Isinfi; (69) 
V2 

0; (70) 
-V2cotn(d u A v - d v A u )) (71) 
Introducing new functions: 

A 2 = A 2 - H cosfl x 3 = A 2 - 6A 2 ; (72) 
lo , 

A v = A v + — H sin fi x 6 = A v - SA V , (73) 
v2 

A 3 =A 3 , (74) 
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where: 

o , 1 o , 

SA 2 =H cosfi x 3 ; 5A V = -= H sin 17 x 3 ; 5A 3 = 0, (75) 

v2 

let us reduce the relations (|68p and to the form similar to (170)) : 

d 3 i 2 - d 2 A 3 = 0; (76) 

9„A 3 - d 3 A v = 0. (77) 

Let us notice that the renormalization of the component of the vector potential (|72[) . 
(|73p keeps the relation ((70)) invariable. But then it is possible to write down the system 
of equations ((701), ([75]). ([TTj) as: 

- = 0; (<r,* = «,2,3) (78) 

and to consider it as equations on a three-dimensional hypersurface V 3 = {v, x 2 ,x }. 
As it is known, the unique solution of equations (|T8[) on V 3 is a gradient function: 

A CT = 9 CT $, (o- = »,2,3), (79) 

where $ = 4>(u, w, x 2 , a; 3 ) is an arbitrary scalar function. The value of the potential 
function corresponding to the initial conditions (|56p is: 

$(u < 0) = x 3 H ( —j={ v — u) sin Q — x 2 cos fl^j . (80) 

Thus 

A rT = d^ + 5A rT . (81) 

As it is known (see, for example, |10|). it is possible to impose one gauge 
condition on 4 components of a vector potential. We choose this condition in the 
form corresponding to the initial conditions (|56[) : 

A u = 0. (82) 
Then for the nonconserved components of the Maxwell tensor F ur7 is valid: 

F w = 8 wr $; (<t = v,2,3). (83) 
But then condition (|7ip can be written down in the form: 

d u (A 2 + %/2cot QA V ) = 0. (84) 
Integrating ([84]) with the initial conditions ([56]) . we obtain: 

A 2 + \/2cotO^ = 0. (85) 
Taking into account the identity: 

5A 2 +\/2cotn<L4„ = 0, (86) 
we obtain the linear differential equation from (f85|) : 



a 2 $ + \/2cot«a u $ = 0. 

Integrating it, we obtain: 

$ = $(wv / 2sinfi-a; 2 cos^,u,x 3 ), (87) 

where $ is an arbitrary function of its arguments. Using now the initial condition (|80p 
we obtain finally: 

$ = x 3 H (~~^( v ~ sinfi- x 2 cosO 
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where tp(u) is an arbitrary function of the retarded time, satisfying the initial 
condition: 

il>(u <0) = u. (89) 
Thus, the final expression for components of the vector potential becomes: 

A-2 = A v = A u = 0; A 3 =H (^( v ~ ^( u )) sinfi - x2 cosQ^j . (90) 
The components of the Maxwell tensor relative to the potential are equal to: 
F vu = 0; F 2u = 0; F 3u = -j= H if/ sin il; 

F 2v = 0; F 3v = — -= H sin O; F 23 = - H cos SI (91) 
V2 

and are defined only by one unknown function ip(u). For the components of the dual 
Maxwell tensor we get: 

F vu = -4 H cosfl; F 2u = jL- H smtt; F 3u = 0; 

L z \/2L 2 

F 2v = H ^'sinfi; F 3v = 0; i^ 23 = 0. (92) 

V2L 2 

3.5. Accompanying intensities and the frozen- in condition 

According to (|17p . we define the components of the vector of accompanying intensity 

of the electric field, Ei, as: 

lo 10, „ 

E V = -—=H sinO v s ; E u — — = H ijj sinQ, v 6 ; 
v 2 v 2 

£ 2 =H cosfi w 3 ; E 3 = —= H sini}(v v — i/j'v u )— H cosil v 2 . (93) 

v2 

Thus, the frozen-in condition of magnetic field in plasma © reduces to two equalities: 

v 3 =0; -^=(v v tp' - v u )sinQ + v 2 cosQ = 0. (94) 
v2 

As a result, covariant components of the Maxwell tensor, contravariant components 
of the dual Maxwell tensor and contravariant components of the velocity vector are 
defined by the expressions, obtained in [T], but now they are already defined for a 
more general metric of a gravitational wave and at weaker assumptions. In the quoted 
paper, in particular, to obtain the explicit form of the Maxwell tensor components 
and the velocity vector components, the analysis of the drift current components was 
carried out using the conservation law of this current. As it was shown above, for 
achievement of this purpose three assumptions are sufficient: 

1. inheritance of the space symmetry by the energy- momentum tensor of the electro- 
magnetic field and by the energy-momentum tensor of the plasma separately; 

2. the equality to zero of the second invariant of the Maxwell tensor; 

3. the frozen-in condition of magnetic field in plasma. 



Thus, the analysis of the first group of Maxwell equations and initial conditions 
is sufficient. 
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Calculating further the covariant components of the dual Maxwell tensor, subject 
to (JUS), we get: 

oo o 

Fuv= yrrcosft; Fu2— —^e 2 ^ cosh 27?// sin ft; F V 2— ^e 2 * 3 cosh 2"f sin ft; 

L v2 v 2 

o o 

F„ 3 = -p sinh 2 7 sin ft; F«3= 7= sinh 2 7 V>' sin ft; F 2 3 = 0. (95) 

v 2 v 2 

Covariant components of the vector of magnetic field intensity relative to the 

Maxwell tensor (|9"2"|) are equal to: 



H v = — \v v cos ft H — j^v 2 sin 17^) ; (96) 




L 2 ■ ^ 



H u = — [vu cos ft f u V sin ft ; (97) 

L V V2 / 

ff 2 = — 7= H cosh2 7 e 2,3 sinft(u t ,^' + v u ); (98) 
v2 
1 

H 3 = —= H sinh 2 7 sin Sl(v v ip' + v u ). (99) 
v2 

It is thus easy to show on the basis of formula (|9"2"|) : 

H 3 =F t3 v i = 0, (100) 

i.e., the third contravariant coordinate of the vector of the magnetic field intensity, as 
well as the vector of dynamic velocity of the plasma, is equal to zero. Also it is easy to 
be convinced of orthogonality of the velocity vector and the magnetic field intensity 
(HI: 



H t v l EE 0. (101) 

The square of the magnetic field intensity, i.e., a scalar H 2 , most easier to calculate by 
means of the relation (|2T|) , using the explicit form of contravariant (|92|) and covariant 
(|95p components of the dual Maxwell tensor: 

H 2 

H 2 = ^- T (L 2 ip' cosh 2 7 e 2/3 sin 2 ft + cos 2 ft) . (102) 

The frozen-in conditions of magnetic field in plasma establishes the connection 
between nonzero contravariant components of the velocity vector v 2 , v° = v u , v u — v v . 
Besides, there is still the normalization relation of velocity vector ([T]). Therefore, 
the only one independent coordinate of the velocity vector remains, and the 
electromagnetic field is defined by only one unknown function of the retarded time, 
ip(u). Using (|96"l) - (|102p the normalization relation of velocity vector can be written in 
the equivalent form: 



v v cos ft + vo — f= sin ft 
V2 



If 2 v , l4 _ sin^ft l2 cogh ^ (iQ3) 




H 2 
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3.6. Drift current 

Let us calculate components of a drift current, using Maxwell equations (|10[) . 
considering the dependence of Maxwell tensor components only on the retarded time 

1 



Then: 



■4 = ~^d u (L 2 F™). (104) 
J£ = J& = 0; (105) 

Ji-A^^.y, doe) 



Ji = - 2y/ 2 nL 2 ( sinh2 ^ • V + cosh2 ^ • ^) • ( 107 ) 

Calculating scalar product of the vector of the magnetic field intensity and the vector 
of drift current density, using ([55]). (|M|) . (ITU5|) . (JTOSJ), ([TOT]) , we get: 

(Jd*,H) = ^j(«„0' +« u )(y - ^sinh4 7 ). (108) 

Thus, the presence of the second polarization of a gravitational wave leads to 
violation of the orthogonality of the vectors of drift current density and magnetic 
field intensitjEl 

Using expression (TTTj) . it is possible to show that the equality (|105|) is carried out 
only in the case of transverse propagation of the PGW (Q = ir/2). 

3. 7. Integrals of the motion 

Because of existence of the motions (|4"6")l . Killing equations are satisfied: 

i £ fcl i = 0> (a = 373). (109) 

(a) (a) 

Therefore conservation laws of the total EMT in a field of PGW after consistently 
transvection with all Killing's vectors (|46[) can be written down in the form: 

1 ' " ' - irpk 



(0*V=S ri?) = 0; (a = 1,3). (110) 



V-9 (a) 

Taking into account the fact that EMT components can depend only on the retarded 
time, we obtain following integrals [3J: 

L 2 £ l T m =C a = Const; (a = M) . (Ill) 

(a) 

In this paper we consider only the case of transverse propagation of the PGW (fi = 
7r/2). Then, substituting expressions for the EMT of the plasma and electromagnetic 
field in the integrals (II 1 1 [) . using relations ff99]) - (| 103[) and also initial conditions (|43|) . 
(|44l) . we lead integrals of the motion to the form: 

2L 2 ( £ +p||)^-(p|| -p ± )^cosh2 7 e 2 ^(e+P)A( U ); (112) 

+ We remind that in case of monopolarized gravitational waves these spacelike vectors are mutually 
orthogonal pp. 
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L 2 (e+ Pl] )v v v 2 =0; (113) 

L 2 (e+ Pll )v v v 3 = 0, (114) 

where: 

p=p± ; (115) 

and so-called the governing function of GMSW is introduced: 

A(u) = 1 - a 2 (cosh2 7 e 2 ' 3 - 1) , (116) 
with dimensionless parameter a 2 : 

a2 = 77^v 11171 

4tt(£ + P) 

Solving (|112l) with respect to v v we obtain expressions for coordinates of the velocity 
vector as functions of the scalars : e, py, p±, ip' and the explicit functions of the 
retarded time: 

2 (g+jj) . . . (P|| -p±) g 2 cosh2 7 e^ 
V " = 2L 2 (e+ P ||) A(M) + TTT^^^Z 2 — 5 (U8) 

From ([ITll we get: 

v 2 =0- (119) 

We obtain the coordinate v u from a normalization relation of velocity vector, using 

cud, en : 

«« = ^" - (120) 

2v v 

and from the frozen-in condition (|M|) we get the value of a derivative of potential ij;': 
1 

using it, the scalar H 2 is defined from relation (|102p as: 

2 _ jf 2 cosh 2 7 e 2 ^ 
H - ^ 2v 2 ■ (122) 
Let us notice that in the case of an isotropic plasma (pj_ = p\\ = p) the expression 
(|118p becomes: 

»« = J(FT15 A <»» i (123) 

From RMHD system of equations it is possible to obtain a following differential 
equation in the PGW metric: 

L 2 e'v v + {e+p\\){L 2 v v )' + ^L 2 {p\\ - Pl _)v v {\nH 2 )' = 0. (124) 

Finally, the equation (|124[) is the differential equation on 3 unknown scalar functions: 
£, P\\ and p±. Such underdefiniteness is a known consequence of the incompleteness 
of hydrodynamic description of a plasma. To solve this equation it is necessary to 
impose two additional connections between functions e, py, pj_, i.e., an equation of 
state: 

P||=/(e)l P± = <?(£)• (125) 



^' = — 2 , (121) 
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4. Barotropic equation of state 

General formulas 

Let us consider a barotropic state of an anisotropic plasma, when the connections 
(|125|) are linear: 

p\\=k\\e\ p± = k±s , (126) 

The equation (|124[) is easy to integrate at the connections (|126p , and we get one more 
integral: 

e(V2L 2 v v ) ( - 1+k ^H {k ^ k± ' ) =eH {k ^~ k±) . (127) 

Thus, formally the problem is solved, as it is reduced to the solution of the algebraic 
equations system which, however, is still too difficult to solve and analyse. The solution 
is essentially defined by two dimensionless parameters: k± and k\\ . Further we consider 
the special cases of these parameters. 

4-2. Transverse propagation of the PGW 

In the case of a barotropic equation of state at the connections (I126P substitution of 
(fT2"2"]) in (fTT5|) leads to result: 

vl = \^{u). (128) 

Substituting (I122[) , (|128[) in (|127p , we obtain the closed equation relative to the variable 
e, solving which, we get definitively: 



o 

e =e 



A l + fc_ Li 2(l+/c || )( cosh27e 2/3> ) fc || -fc_ L . ( 12 g) 



-g± 



1 

V2 



AL (fc||+^)( cogn27e 2/3)- 



9 A 







AL( 1+fe »)(cosh2 7 e 2 ^- — 



(130) 
(131) 



where 



£1=^-^-6 [1,2]. (132) 
In particular, for ultrarelativistic plasma with zero parallel pressure: 

fe||->0; k ± ^^ (133) 
we obtain from (|T29l) - (fT32|) : 

v v = ^LA 2 (cosh2 7 e 2 ^)- 1/2 ; (134) 
V 2 

e=eL- 4 A- 3 (cosh2 7 e 2/3 ); i? =1 L- 2 A- 2 (cosh 2 7 e 2/3 ) . (135) 
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5. Energy balance equation 

In [1] was shown that the singular state, which exists in a magnetized plasma under 
the condition 2/3oa 2 > 1 on the hypersurface: 

A(u*)=0, (136) 

is removed using the back effect of the magnetoactive plasma on the GW. That leads to 
the efficient absorption of GW energy by the plasma and restriction on the amplitude 
of the GW. An exact solution of the PGW energy transformation to the energy of 
the shock wave is possible only on the basis of the self-consistent system of Einstein's 
equations and magnetohydrodynamics equations. However, qualitative analysis of this 
situation can be carried out using a simple model of energy balance proposed in [2] . 
The energy flow of the magnetoactive plasma is directed along the direction of the 
PGW propagation, i.e., along the axis Ox 1 . Let /3*(u) and 7*(u) are the vacuum 
amplitudes of the PGW. In WKB-approximation: 

8tt£ < uj 2 , (137) 

where ui is the characteristic PGW frequency and e is the matter energy density, all 
the functions still depend only on the retarded time (see [12])- Thus, (3(u) and ^(u) 
are the amplitudes of the PGW subject to absorption in plasmas; is the total 
energy-momentum tensor of the plasma and the electromagnetic field ([2} . 



5.1. Integral law of energy conservation 

Ref. pQ suggested a semiquantitative solution of this problem on the basis of a simple 
model of energy balance. Due to its extreme importance, we do not restrict ourselves to 
PQ and return to a more complete study of the problem of energy transmission from 
a GW to magnetoactive plasma. However, instead of solving Einstein's equations, 
we make use of their consequence, the conservation law of the total momentum of 
the system "plasma + gravitational waves" . Clearly, this model is only approximate 
and cannot replace a rigorous solution of Einstein's equations. According to |10) . an 
arbitrary gravitational field provides the conservation of the system's momentum: 

P l = \J (-9)(T 2i + T li )dV, (138) 

g ., 

where T is the energy-momentum pscudotcnsor of the gravitational field and the 
integration covers the whole 3-dimensional space. Let us take into account that 
the above solution is plane-symmetric and depends on the retarded time u only. 
Consequently the integration over the "plane" (x 2 ,x 3 ) in (1138j) reduce to simply 
multiplying by an infinite 2-dimensional area. Dividing both sides of (|138|) by this 
area and bearing in mind that with Q = ir/2 among the 3-dimensional flow only P 1 is 
nonzero, we obtain the conservation law of the surface density of the momentum P^: 

+oo 

= - c J {-g){T u +T 1A )dx (= Const). (139) 

— oo 

Let the right semispace x > be filled with magnetoactive plasma and the left one 
x < with matter which does not interact with a weak GW. Let further the whole 
gravitational momentum be concentrated in the interval u € [0, Uf] where tf = 
is the gravitational pulse duration. Since the integral in Eq. (1139j) is conserved all the 



Exact solutions of relativistic magnetohydrodynamics equations 17 

time, let us consider it at to < 0, when the GW has not yet reached the magnetoactive 
plasma, and — tt > t > 0), when the GW has reached the plasma. Taking into account 
that the vacuum solution depends only on the retarded time, we get for the integral 
in Eq. JBSJ: 

u f ts/2 u f 

T l A du = J (T 14 + T u )du + J T l A du, (140) 

t/V2 

where T \ A =T 14 (/3*(w),7*(w)); T 14 =T 14 (/3(w), 7 (u)). Transferring one of the 
integrals to the left-hand side of Eq. (|140l) . we obtain the relation: 



Tt?du = J (T + T )du, (141) 
o o 

where the variable u = t/^/2 > can now take any positive values. 

A similar law may be written for the plasma total energy; in this case instead of 
Eq. (|14ip we obtain: 



J Tfdu = |(T 44 -f +T 14 )d U , 



o o 

where E o is the total energy density of the unperturbed plasma. 
5.2. Local analysis of the conservation law 

Since the relation (|141l) must be valid at any values of the variable u, the corresponding 
local relation should be satisfied: 

T 41 (/3,7)+T 41 (/3,7)=T 41 (/3*,7*), (142) 

g 

where T is the energy flow of a weak GW in the direction Ox 1 (see |10|): 



9.41 1 



T 

16tt 



1 



1 



hg + - A (K 2 - h' 33 f =—[( 1 'f + (P'f\. (143) 



4 s 

The prime denotes differentiation with respect to s. At substituting (|143[) into (|142p 
and changing variables to v, u, we obtain: 



2tt 



rpVV rjiUU 



In the case of transversal PGW propagation and at a barotropic equation of state of 
an anisotropic plasma we obtain: 

r~ - T wu =(Jg- vlj (e(l + fc ± ) + |J) . (145) 

Further, using the solutions of magnetohydrodynamics for a barotropic equation of 
state of plasma (|129[) . (|130[) . (|131[) and dimensionless parameter a 2 (|117p . we rewrite 
the energy balance equation (11441) as: 

jj2 / 4 4(fc||+fcj_) 2(fcy -fejj \ / 1 \ 

— (A (cosh2 7 e 2/3 ) ^ - 1 j ( — + 1 j + 

(V) 2 + = (Y,) a + (#) a . (146) 
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Let us expand the expression in brackets by the smallness of the PGW amplitudes 
(|4ip but hold the term with A -1 , since the parameter a 2 in a strongly magnetized 
plasmas can be so large that the condition 2a 2 /3 > 1 is satisfied. Then energy balance 
equation takes the form: 

^ (A- 49i - 1) + l) + (V) 2 + (/3') 2 = «f + m 2 - (147) 

Since in linear approximation by the smallness of the amplitudes (3 and 7 the governing 
function (| 1 1 6[) does not depend on the function 7(1*): 

A(tt) = l-2a 2 /? + 0(/3 2 , 7 2 ), (148) 

and the functions f3(u), "f(u) are arbitrary and functionally independent, then up to 
/3 2 ,7 2 , the relation (|147[) should be decompose into two independent parts: 

2H 2 g ± (l + a 2 )(3 + (P') 2 = m 2 , (149) 
(l') 2 = «f- (150) 
Here, according to the meaning of local energy balance equation, we consider short 
gravitational waves (|137|) . so we can neglect the squares of the PGW amplitudes in 
comparison with the squares of their derivatives with respect to the retarded time. 
Thus, according to (|150l) : 

7*(u)= 7 (u), (151) 

i.e., in the linear approximation a weak gravitational waves with polarization e x does 
not interact with a magnetized plasma. This coincides with the conclusion of the 
paper [13]. 

Thus, in the linear approximation the PGW with e x polarization passes through 
a magnetoactive plasma without absorption, and the energy balance equation takes 
the form obtained in [6] . Further conclusions are similar to the case of propagation of 
the PGW with only one polarization e+. 

If a 2 3> 1 the Eq. (| 14T[) can be written in the form (see also [2]): 

eV{q) +q 2 = q 2 , (152) 

where q — (3//3q, the dot denotes differentiation with respect to dimensionless time 
variable s: 

s = V2uu, (153) 
(to - the PGW frequency), V(q) - potential function which in a weak PGW becomes: 

V(q) = A- 4 ^(q)-l, (154) 
where £ 2 is so-called the first parameter of GMSW [2J: 

r2 



Let us introduce the new dimensionless parameter: 

T = 2a 2 /3 (156) 
- (the second GMSW parameter) and rewrite (I148[) in a weak PGW as: 

A(q(s)) = 1 - 2a 2 f3 Q q(s) = 1 - Tq(s). (157) 
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It leads from (|T57|) : 

4—^. (158) 

To analyze the system behavior, let us suppose that the moment s = corresponds 
to the front edge of a GW, while: 

ft « /3 (1 - cos(s)) =S> w 1 - cos(s). (159) 

According to (|157|) - (I159|) the system starts with negative value of the governing 
function derivative and with function value equal to 1: 

A(s) ps — Tsins rj — Ts; 

2 (s -> +0). (160) 
A(s) « 1 - T(l - coss) « 1 - Ty; 

The energy balance equation (I152[) according to (|154|) . (I158[) . (I159P becomes: 

A 2 + £ 2 T 2 [a~ 4 ^ - lj = T 2 sin 2 (s). (161) 

The minimum value of the governing function at s — 7r/2 is equal to: 

1 \-7x 



where: 



A mm = (— + !] , (162) 



1 l-fcj. 1 1 „„„s 

^ = 4^1 = — ^8^1- (163) 



The maximum accessible density of a magnetic energy is 

h 2 \ h 2 I r 



and it does not depend on a plasma equation of state. Also the plasma velocity in a 
GMSW does not depend on an equation of state. And the maximum plasma energy 
density without magnetic field depends on the exponent of plasma anisotropy: 

e max =e I 1 + (165) 
It is maximum for ultrarelativistic plasma with zero parallel pressure. 



6. Conclusion 



Thus, the generalization of the results of [T] -[3] in the case of gravitational wave 
with two polarizations has been obtained and has been shown that in the linear 
approximation the polarization e x does not interact with a magnetized plasma. This 
fact is a justification for applicability of the previously obtained results for the case of 
arbitrarily polarized gravitational wave. 
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